Intermittency of velocity time increments in turbulence. 
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We analyze the statistics of turbulent velocity fluctuations in the time domain. Three cases are 
computed numerically and compared: (i) the time traces of Lagrangian fluid particles in a (3D) 
turbulent flow (referred to as the dynamic case); (ii) the time evolution of tracers advected by a 
frozen turbulent field (the static case), and (iii) the evolution in time of the velocity recorded at 
a fixed location in an evolving Eulerian velocity field, as it would be measured by a local probe 
(referred to as the virtual probe case). We observe that the static case and the virtual probe cases 
share many properties with Eulerian velocity statistics. The dynamic (Lagrangian) case is clearly 
different; it bears the signature of the global dynamics of the flow. 



One of the distinctive feature of turbulence is the de- 
velopment of extremely high fluctuation level at small 
scales. The probability density functions (PDFs) of ve- 
locity increments are stretched at small scales, while they 
are almost Gaussian at large scale where energy is fed 
into the flow Q. This evolution is traditionally referred 
to as "intermittency" in turbulence studies. Numerous 
studies have been devoted to the study of intermittency 
in the spatial domain, analyzing velocity differences be- 
tween two points separated by a variable distance r. For 
instance, it is now well established that when the velocity 
increments are computed along the distance r, the longi- 
tudinal velocity increments are skewed and the structure 
functions have universal relative scaling exponents in the 
inertial range 0, 0, 0] , related to the graininess of the 
dissipation. When the increments are related to changes 
in the velocity component perpendicular to the distance 
r, the corresponding transverse structure functions dis- 
play a different scaling, related to the spatial distribution 
of vorticity 0- In contrast, there has been much fewer 
studies of intermittency in the time domain, i.e., related 
to changes in time of the velocity field. Two cases are of 
particular importance: The Lagrangian one, which per- 
tains to the fluctuations in time of the velocity of marked 
fluid particles, and the Eulerian one, where one consid- 
ers the variations in time of the velocity measured at 
a fixed location in the flow. We consider of course the 
case where turbulence develops in the absence of a mean 
flow, otherwise Taylor's hypothesis trivially reduces the 
second case to a spatial measurement Eulerian time 
fluctuations are different because one expects, after Ten- 
nekes original suggestion Q , that sweeping (the random 
advection by large scale motion) plays an important role. 
In practice, Eulerian intermittency in time is relevant for 
stationary bodies exposed to turbulent flow conditions 
(atmospheric ones, for instance). Lagrangian intermit- 
tency, on the other hand, has strong implications for pro- 
cesses such as mixing , combustion p , or cloud forma- 
tion [ToT | . Lagrangian data has recently been made avail- 
able in numerical simulations [ill . fl2T | and experimental 
measurements ^1 Q, ^| . One rather unexpected fea- 



ture is the observation of long-time correlations in the 
Lagrangian dynamics ^(| . It has been incorporated in re- 
cent stochastic models of Lagrangian acceleration [TtIIisI 

- see also for a recent review. However, the full 
modelization of Lagrangian velocity fluctuations is still in 
progress. Lagrangian statistics are often computed from 
the advection of particles in a frozen Eulerian field (to 
minimize computing overhead). We have thus decided 
to compare this pseudo-Lagrangian statistics to that of 
pure Lagrangian and Eulerian time fluctuations. 

We study the problem numerically. We first compute 
the velocity changes of fluid particles that are advected 
by a frozen 3D Eulerian velocity field; that is, we consider 
a single snapshot of a converged turbulent flow, and use 
it to advect fluid particles in this frozen Eulerian field 

- we call this the static case. Then we compute the 
velocity variations of true Lagrangian particles, a situa- 
tion in which the Eulerian flow is also evolved in time 
according to the Navier-Stokes equations — this case is 
called the dynamic case. Finally, we record the time evo- 
lution of the Eulerian velocity at fixed locations of the 
computation domain, as it would be measured by virtual 
velocity probes. We then perform a comparative study 
of the intermittency characteristics of the three velocity 
time signals. We show that the statistics of time veloc- 
ity increments depends on the situation considered: the 
static case and the virtual probes display intermittency 
features that are reminiscent of Eulerian velocity statis- 
tics; the former has multifractal spectra identical to the 
ones measured for 3D Eulerian turbulence jl^l , while the 
later coincides with traditional longitudinal Eulerian ve- 
locity increments statistics Q . The dynamic case, which 
is a true Lagrangian measurement, displays significantly 
more intermittent features. 

The Navier-Stokes equations are integrated in a 256 3 
cubic domain of size 2-k by a parallel distributed mem- 
ory pseudo-spectral solver, using a second-order (in time) 
leap-frog scheme. The large-scale kinetic-energy forcing 
is adjusted at each time step by scaling the amplitudes 
of modes 1.5 < k < 2.5 uniformly (phases are left to 
fluctuate freely), so as to compensate exactly the losses 
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due to eddy dissipation in the kinetic-energy budget [2l| . 
The Reynolds number based on the Taylor microscale is 
R\ = 140. The velocity root mean square is 0.1214 m/s, 
the mean dissipation rate is e = 0.0011m 2 /s 3 and the 
kinematic viscosity is v = 1.5 10 -4 m 2 /s. The particles 
trajectories are resolved in time by a second-order Rungc- 
Kutta scheme, and interpolated using cubic spline func- 
tions. In the dynamic case, a set of 10,000 particles, uni- 
formly distributed in the cube at initial time, are followed 
for a duration of approximatively 7Tj,. Here, as in |23l |. 
Tl is defined for each signal as the time scale above which 
the velocity statistics is Gaussian (the second order cu- 
mulant has reached the Gaussian value 7r 2 /8). For the 
static case, 10,000 trajectories have been integrated over 
approximatively 3Tl- Finally, 32,768 virtual probes have 
been used to get the time variation of the Eulerian field 
at a fixed location; in this case the records are 10 Tl long. 
We label vo,i(f) the Lagrangian velocity of one compo- 
nent (in cartesian coordinates) of particle number i in the 
3D Eulerian time-evolving flow (dynamic case), vs,i(t) 
one component of the velocity advected by the static Eu- 
lerian flow (static case) and VT,i(t) the time evolution of 
one component of an Eulerian velocity probe. 

In Fig. n wc show the power spectral densities 
(\v s (uj)\ 2 ), (\v D (u})\ 2 ) and (\v t (uj)\ 2 ) vs. ujT l in a loga- 
rithmic representation, where v means Fourier transform 
of v, for the static, the dynamic cases (averaged over 
the number of tracked particles) and for the virtual fixed 
probes (averaged over the number of probes). One ob- 
serves for all scaling behavior on a small range of 
scales. We show in the inset the values of the correspond- 
ing power law exponent determined from the local log- 
arithmic slope of the spectra. The dynamic Lagrangian 
velocity spectrum has an u>~ 2 inertial range spectrum, 
as expected from Kolmogorov similiraty arguments |24j , 
and in agreement with previous experimental and numer- 
ical observations [l5| . The spectrum of the time variation 
of the Eulerian velocity field (as recorded by the virtual 
probes) shows a clear —5/3 scaling exponent, character- 
istic of Eulerian data Q] . This is in good agreement with 
Tennekes sweeping argument 0, H^): the characteristic 
velocity fluctuations at scale t is the standard deviation 
of the flow velocity v lms = \/ (|fs| 2 ) rather than the Kol- 
mogorov one vi = (et) 1 ' 3 . As a result, a time-scale in- 
crement r corresponds to a length I = w rms T, so that the 
scaling of the increments is in fact Eulerian. In addi- 
tion, this effect produces a larger inertial range because 

3/2 

the dissipative time scale is TlR x for the Eulerian 
time data, while it is TlR^ 1 for the dynamic Lagrangian 
data. Fig. ^ then shows that the data in the static case 
also follow a —5/3 scaling law, closer to the Eulerian be- 
havior than to the Lagrangian one. This result has also 
been found in a similar study using Kinematic Simula- 
tions of turbulence |2^| , a situation in which the Eulerian 
velocity field is monofractal, i.e. not intermittent. Let 
us mention that since inertial range statistics are likely 
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FIG. 1: Power spectral density E(ojTl) of one component of 
the velocity of tracer particles in the static (□) and dynamic 
(o) cases, and for the virtual Eulerian probes (o). Spectra 
have been shifted vertically for clarity. The local slopes of 
these spectra are plotted in the inset. 

to be independent on the Reynolds number, we consider 
the behavior of our estimators in the inertial range (i.e, 
power spectra in Fig. 1 and cumulants of magnitude in 
Fig. 2) as characteristics of fully developed turbulence 
- despite the moderate Reynolds number value of our 
DNS (i? A ^ 140). 

We now seek to quantify the intcrmittency features. 
This is usually done via the analysis of the scaling behav- 
ior of velocity structure functions S(p,r) = (\S T v\ p ) = 
(\v(t + t) — v(t)\ p ) ~ r^p, where the average is computed 
over all accessible times t and over all recorded time se- 
ries. Note that we use the absolute value of velocity 
increments in the definition of the structure functions in 
Eq. since the statistics that are studied are symmetric 
under the transformation 5 T v — > — 8 T v : velocity incre- 
ment statistics are not skewed. However, as advocated 
in [23 for Eulerian velocity data analysis, the magnitude 
cumulant analysis provides a more reliable alternative to 
the structure function method. The relationship between 
the moments of \S T v\ and the cumulants C„(r) of In \5 T v\ 
reads 

(00 n \ 

E C «(^J • (!) 

Previous studies 0, 0, H^] have shown that intermit- 
tency is suitably described for small p using a log-normal 
statistical framework corresponding to a quadratic ( p 
spectrum : 

P 2 

C P = C1P-C2— , (2) 
where the parameters c\ and c-i can be extracted from the 
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FIG. 2: Cumulants of velocity magnitude for the static (□) 
and dynamic (o) cases, and for the virtual Eulerian probes 
(o). (a) Cf D ' T (Y) vs. ln(r/Ti); the solid and dashed 
lines correspond to the slopes cf' T = 1/3 and of = 1/2 
respectively; for the sake of clarity, we have substracted 
{\u\5 Tl v\). (b) C* 2 s,d ' t (t) vs. \xi(t/T l ); the solid, dashed- 
dotted and dashed lines correspond to the slopes cf = 0.046, 
cf = 0.03 and cf = 0.08 respectively; we have subtracted 
Var[ln \8t l v |] = 7r 2 /8 and also shifted the upper curves by 0.1 
and 0.2 for clarity. 



time scale behavior of the first two cumulants C\ (r) and 
Ci(t). In Eulerian context, the analysis of longitudinal 
velocity increments in the incrtial range has shown that 
QfW = — C 2^ n (£/L), where £ is a spatial scale and L 
the decorrelation length, with a universal intermittency 
coefficient cf w 0.025 HH^. 

We report in Fig.[21the results of a comparative analy- 
sis of the cumulants C 1 ' D ' T (t) and C 2 ' D ' T (t). The pro- 
file of Cf (r) as a function of 1ii(t/Tl) is significantly 
curved, a feature that has also been observed on exper- 
imental data of Lagrangian velocity structure functions 
[l5j . This departure from scaling is a signature of (i) 
the pollution of the incrtial range by d issip ative (finite 
Reynolds number) effects as studied in |23j and (ii) the 



non universal and/or anisotropic behavior of the turbu- 
lent flow at scales of the order of the integral time scale 
~ Tl. In Fig. 2(a), we have indicated by a dashed 
line the scaling behavior Cf (r) = cfln(r/Ti), with 
cf = (h) ~ 1/2 corresponding to a (ujTl)~ 2 scaling 
region in the power spectrum (we recall that h ~ 1/2 
corresponds to the most probable velocity scaling expo- 
nent in a multifractal analysis of Lagrangian intermit- 
tency 0,H^|). In contrast, the first order cumulant for 
the static case and for the Eulerian velocity time varia- 
tions is better represented by C x ' (r) = c 1 ' 1ii(t/Xl), 
with c 1 ' T = (h) f=a 1/3, as indicated by the continuous 
line in Fig. |2{a). Again, this value is consistent with the 
(luTl)~ 5 / 3 power spectrum behavior observed in Fig. ^ 
As it is also the case for Eulerian fields 0,H3, the second 
order cumulants — Fig. [3(b) — has a logarithmic bchav- 
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ior in the incrtial range, C 2 ' ' (r) = — c 2 ' ' ln(r/Ti). 
In the dynamic case, we get cf — 0.08, in good agree- 
ment with previous experimental data |23j and numerical 
simulations 0. 

The static case and the Eulerian time probes follow the 
same behavior, with cf ~ 0.046 and cf — 0.03. The first 
value is in good agreement with the recent estimate of 
the intermittency coefficient c^' 30 — 0.049 of a numer- 
ical 3D Eulerian velocity field [2jj. Note that c^ D is 
significantly larger than the value cf — 0.025 computed 
for longitudinal velocity increments, to which one should 
compare the estimate of the intermittency coefficient cf 
of the time variations of the Eulerian longitudinal ve- 
locity component. As a technical but important point, 
we stress that the values of C2 reported here are not com- 
puted from a linear regression fit in the Ci (t) curves (due 
to the narrowness of the inertial range) but as a result of 
the mutifractal description of the entire range of scales, 
dissipative domain included, as detailed in |23j. 

We thus conclude that the static and Eulerian inter- 
mittencies in time have identical statistics to respectively 
3D and ID Eulerian fields, while the true dynamic (La- 
grangian) case is clearly different and more intermittent. 
This is confirmed in Fig. [3] where the more familiar struc- 
ture function exponents ( p are compared. Note that in 
order to be able to compare Eulerian and Lagrangian 
data, all exponents are computed using the second order 
structure function as a reference (extended self-similarity 
(ESS) ansatz [i^|) in order (i) to overcome the observed 
bending of the structure functions when plotted versus 
the time scale in a logarithmic representation (as already 
noticed in Fig. |2Ja) for the first cumulant Ci(r)), and 
(ii) to give a clear picture of intermittency effects respon- 
sible for departure from a linear behavior (i.e. £p ^ p/2 
and CjiCp 7^ p/3)- Thus, in this representation, each 
Cj spectrum must be compared to the monofractal Kol- 
mogorov's prediction ( p = p/2. As the size of the statis- 
tical ensemble available is limited, we restrict ourselves 
to moments of order p up to 6. The values obtained 
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FIG. 3: Structure function exponent Cp/& vs - P as com- 
puted for the static case (□), the dynamic (o) case and 
the Eulerian time probes (o). Are also shown for compar- 
ison the ex pon ents obtained for Lagrangian velocity exper- 
iments (•) Hal. 3D Eulerian velocity fluctuations obtained 
by DNS (V) 1201 . experimental ID Eulerian longitudinal ve- 
locity increments (A) Q and passive scalar increments (*) 
l2iill . The solid, dashed-dotted and dashed lines correspond 
to the quadratic £ p spectra (Eq. (2)) with the parameter val- 
ues: (cf = 1/3; cf = 0.046), (cf = 0.3; tf = 0.03), and 
(cf = 1/2; cf = 0.085). 



are in agreement with a parabolic spectrum (Eq. (0)) 
when fixing the parameters c\ and C2 to the values pre- 
viously estimated from the magnitude cumulants and 
auto-correlation functions. In Fig. [3] are also shown for 
comparison the £ p spectra obtained for experimental La- 
grangian velocity measurements, experimental Eulerian 
longitudinal increments, and full 3D numerical Eulerian 
velocity multifractal analysis. Once again, the static and 
Eulerian time behavior are identical respectively to that 
of the 3D numerical Eulerian velocity |2(| and to the 
traditional ID longitudinal velocity increments Q. But 
they are both less intermittent than observed for the dy- 
namic case, i.e. for the Lagrangian velocity field. A de- 
tailed account of the relationship between Eulerian and 
Lagrangian intcrmittencies, in the framework originally 
proposed by M. Borgas j2^| has been previously discussed 
in [23j. Finally, we note that another quite noteworthy 
feature in Fig. [3] is that the exponents for the passive 
scalar increments are identical to that of the Lagrangian 
velocity statistics (within error bars). 

To summarize our findings, we have observed that the 
statistics of particles advected in a frozen Eulerian field 
is, to some extent, similar to that of the time variations 
of the Eulerian velocity at fixed points in space. This 
is an crgodicity property of homogeneous, isotropic tur- 
bulence. The similarity of the static case and the full 
3D Eulerian field could prove useful because local varia- 



tions in time are way easier to measure than full spatial 
3D flows — although it is absolutely necessary that the 
mean flow be truly absent, otherwise local time varia- 
tions relate to spatial profiles, using Taylor's hypothe- 
sis ||1 01- Finally, the intermittency measured in the 
dynamics case is different, showing that true Lagrangian 
particles are sensitive to the global time evolution of the 
flow. One eventually expects that the large scale dynam- 
ics is even more crucial in understanding mixing effects 
in real non-homogeneous flows. 
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